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ABSTRACT 

We present the results of dust scattering simulations carried out for the Orion Eridanus Superbubble region 
by comparing them with observations made in the far-ultraviolet. The albedo and the phase function asymmetry 
factor (g-factor) of interstellar grains were estimated, as were the distance and thickness of the dust layers. The 
results are as follows: 0.43^g for the albedo and 0.45^q 2 for the g-factor, in good agreement with previous 
determinations and theoretical predictions. The distance of the assumed single dust layer, modeled for the 
Orion Molecular Cloud Complex, was estimated to be ~'110 pc and the thickness ranged from ^130 at the 
core to ^50 pc at the boundary for the region of the present interest, implying that the dust cloud is located in 
front of the superbubble. The simulation result also indicates that a thin (^10 pc) dust shell surrounds the inner 
X-ray cavities of hot gas at a distance of ~70-90 pc. 

Subject headings: ISM: bubbles — ISM: dust, extinction — ISM: individual (Orion-Eridanus) — ISM: struc- 
ture — ultraviolet: ISM 



1. INTRODUCTION 

As stellar ultraviolet (UV) radiations scattered by interstel- 
lar dust grains are regard ed as the most d ominant source of 
the diffuse Galactic light dSeon et al.ll201 Ih . the optical prop- 
erties of the dust grains, generally characterized by the albedo 
and the phase function asymmetry factor g = < cos 9 >, can 
be inferred from the observations of the diffuse Galactic light. 
The observations can also be compared with model calcula- 
tions. For example, Draine (2003) estimated the albedo and 
g-factor for the carbonaceous-silicate grains, and obtained 
^0.40 and '^0.65 for the albedo and the ^-factor at A ^ 
1600 A, respectively. While the values of the albedo and 
g-factor are clearly dependent on the models of dust grains, 
accurate photometry and modeling of dust clouds are also re- 
quired for reliable estimations of these optical properties from 
observations. Until now, there have been numerous stud- 
ies in this line of research from the early 70's, especially at 
far-ultraviolet (FUV ) and n ear-ultraviolet (NUV) wavelengths 
jWitt & Lillie 1973, 1978; Lillie & mnl9W; Morgan et al. 
1978, 1982; Morgan 1980^ Onaka et a l, 198 4; Witt_^et al. 
1992, 1997; Gordo n eTaTI [l994: Calz ettTetan [T995h . For 
example, the observations made from the Voyager 2 space- 
craft have been used for the wavelengths shorter than 1200 
A ("Witt et al. 1993; Murthv et al. 1993; Burah et al." 2002t 
[ShaHma & Murthy 2004; Shalima et al. 2006^^ Sujatha et al] 
llOOS, 2007). A comparative review of past work can be found 
in iGordon (.200 4). Recently, GALEX observations provided 
diffuse UV background images with high spatial resolution, 
which have been used for the determination of the albedo and 
the g-factor in both the FUV (1350-1750 A) and NUV (1750- 
2850 A) bands ( Sujatha et al. 2009, 2010). Most of these stud- 
ies showed that the interstellar grain causes a strong forward 
scattering in the FUV, with a moderate albedo, which is in 
good agreement with theoretical predictions. 

The Orion-Eridanus Superbubble (henceforth OES) is a 
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large hot bubble with a diameter of ^30° across the sky 
located in the Orion and Eridanus constellations. Its dis- 
tance is known to be ^155 pc to the near-side and ^^586 
pc to the far-side of the bubble, and its near edge is consid- 
ered to interact with the expanding shell of the Local Bub- 
ble dBurrows et al.ll 19961) . The OES has been observed exten- 
sively in various wavelengths, ranging from X-ray to radio. 
Enhanced soft X-ray emission seen at its inner cavities is the 
evidence of the hot gas (T ^ 10'' K) occupying its inner parts 
(Heiles et al. 1999). Two prominent filamentary Ha shells, 
known as arcs A and B, have been obs erved at the bound- 
aries of these hot cavities (iReynolds et al.iil998: Boumis et alJ 
I2OOI . The presence of 21 cm radiation of neutral hydrogen 
(T ~ 10^"^ K) was also noted outside the hot bubble (Brown 
et al. 1995). The dust extinction {E(B-V)) map, made from 
the Schlegel, Finkbeiner and Davis (henceforth SFD) Dust 
Survey (Schlegel et a l. 1998), shows a strong anti-co rrelation 
with the soft X-ray emission image (iJo et alJl201 ll) . Recent 
spectroscopic observations made in the FUV revealed the ex- 
istence of C IV and Si II* emission lines as well as molecular 
fluorescence H2 lines, implying that the OES is truly a multi- 
phase object, consisti ng of hot gases and cold dust that in- 
teract with each oflie r (iKregenow et al.ll2006l: iRyu et alJl2006l 
IMjStlJcTetaniMlh . 

The dust properties of the OES region have not been ex- 
plored well because of the lack of sufficient diffuse FUV ob- 
servations. Ba sed on the data av ailable from the Voyager 
measurement, Murthv et al.l (11993) suggested a lower limit of 
^0.3 for the albedo assuming isotropic scattering, and an up- 
per limit of ^0.8 for the phase function asymmetry factor with 
perfectly reflecting grains. However, it should be noted that 
the Voyager observations provided limited information, and 
only for two regions with a field of view of 0°. 10 x 0°.87 and 
a spectral band of 38 A. In this paper, we report the results 
of our study on the dust scattering properties of the OES re- 
gion. We performed Monte Carlo simulations and compared 
the results with the diffuse emission map made from the re- 
cent FUV imaging spectrograph mission. We obtained the 
average albedo and phase function asymmetry factor, and es- 
timated the distance and thickness of the dust cloud. 
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Fig. 1. — (a) FUV continuum map overplotted with dust contours, in equa- 
torial coordinates. The continuum intensity is given in CU (photons s"' cm 
sr"' A"') and the levels of dust extinction contours are in E(B-V). (b) A scat- 
ter plot for the FUV intensity against dust extinction, obtained from pixel-by- 
pixel comparison of Figure 
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2. OBSERVATIONS AND DATA REDUCTION 

The present study is based on the same data set ob- 
tained via Far-Ultraviolet Imaging Spectrograph (FIMS) as 
was used previously in the studies of iRvu et alj ( 120061 120081) 
and lJo et alJ dloTlh . FIMS is a dual-channel FUV imaging 
spectrograph (S-band 900 A-1150 A, L-band 1350 A-1750 
A) on board the Korean microsatellite STSAT-1, which was 
launched into a 700 km sun-synchronous orbit on September 
27, 2003 and operated for a year and a half. With its large 
image fields of view (S-band: 4°.0x4'.6, L-band: 7°.5x4'.3) 
and moderate spectral (A/AA ~ 550) and angular (^^5') res- 
olutions, the instr ument is optimized fo r the observation of 
diffuse emissions jEdelstein et al.ll2006allbl) . Though the OES 
region was observed in both the S- and L-band, we ana- 
lyzed only the L-band data for the present study since the 
S-band data were strongly contaminated by geocoronal emis- 
sion lines. Bright point sources were removed from the map: 



a total of 210 stars in the TD-1 catalog with 1565 A band 
fluxes of above 10"'^ ergs A"' s~' cm~^ were removed, and the 
remaining bright pixels, not corresponding to the stars listed 
in the TD-1 catalog, were also removed. More informa tion o n 
the data reduction procedures can be found in Jo et al.l(l2011l) . 

While stellar UV radiation scattered by interstellar dust 
grains is generally known to be the most dominant source 
of diffuse FUV emission, there are other contributors such 
as atomic line emissions, extragalactic background, H2 fluo- 
rescence, and two-photon continuum emissions. Hence, it is 
necessary to eliminate these extra contributions in order to ob- 
tain more reliable results for the effect of dust scattering. For 
example, the atomic lines were removed as follows. We took 
only the 1420 A-1630 A portion from the L-band spectrum to 
exclude the bright airglow lines at O I 1356 A and 1641 A as 
well as other atomic/ionic lines such as O IV] 1407 A, He II 
1640 A, CI 1657 A , and AI II 16 71 A, as identified previ- 
ously in Table 1 of J o et al.l (1201 Ih . The remaining spectrum 
still included the Si II* and C IV lines centered at 1533 A and 
1549 A, respectively, which were removed by excluding the 
spectral range from 1525 A to 1560 A. We have constructed a 
diffuse FUV map with the remaining portion of the spectrum. 

The extragalactic diffuse FUV background is known to be 
more or less uniform over the celestial sphere with ^^300 
Continuum Units (CU), a nd is expected to be even more so 
at high galactic latitudes (lBowveij|199Th . Noting that OES 
is located at high galactic latitudes, we have subtracted 300 
CU from the original diffuse emission map. The contribu- 
tion of H2 fluorescent emission was subt racted in the fo llow- 
ing way using the map of Figure 2(b) in lJo et al.l (l20Tl1) . We 
first converted the original H2 map, which was constructed 
from the 1608 A line intensity, to a map of total H2 emission 
for the whole wavelength range in consideration of using the 
H2 fluorescence emission profiles given by the CLOUD co de 
dvan Dishoeck & Blackll^lBlack & van Dishoeckll987h . a 
simulation model for a photodissociation region (PDR). We 
assumed that the intensity ratios among the lines were the 
same for the entire OES region. The resulting emission in- 
tensity was expressed in terms of CU by dividing the total 
intensity with the wavelength range, and subtracted from the 
diffuse emission map. 

It was mentioned in iJo et al.l (1201 lb that two-photon emis- 
sion constitutes about 50%-70% of the FUV continuum in- 
tensity in the lower part of arc A and arc B, based on the 
assumption that 1 Rayleigh of the Ha line corresponds to the 
two-photon (plus free-free and bound-free) con tinuum emis- 
sion o f 60 CU in an ioni zed gas dRevnoldsl 19901) . On the other 
hand, ISeon et all (1201 ll) recently suggested that the contribu- 
tion of two-photon effects to the FUV continuum could be 
lower in the case of cooling ionized gas with 1 Rayleigh of 
the Ha line cor responding to ^26.0 CU. Furthermore, Figure 
6 of IWitt et al] (I2OIOI). which shows the data for the OES re- 
gion studied bv lMadsen et al.l f2006 ). suggests that 10 to 40% 
of the total Ha intensity comes from dust scattering. These 
complexities make it difficult to estimate the effect of two- 
photon emissions. Therefore, we have not removed the two- 
photon continuum emission, which may make the two arcs 
appear more prominently, and will instead discuss its contri- 
bution to the map. The final image was obtained with 1° x 1° 
pixel bins and is shown in Figure ^a), overplotted with dust 
contours. 

When the new map of Figurellla) is compared with the pre- 
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vious diffuse FUV map. Figure 2(a) of lJo etaUdMB . which 
includes the contributions of the ion lines as well as the H2 
fluorescent emission, we can see that overall morphology did 
not change much. The most significant change due to the ex- 
clusion of the ion lines is seen only in the northeastern region 
of the map near the left boundary above S = 0°, where the 
scattered photons of the two bright B-type stars with strong P 
Cygni profiles, HD 31237 and HD 30836 located at (a, 6) - 
(73°.6, 2° A) and (72°.8, 5°.6), respectively, are responsible 
for the enhanced C IV emissions. The H2 fluorescent emis- 
sion varied significantly according to the regions as can be 
seen in Figure 2(b) of iJo et all dlOl iT) while it is -17% on av- 
erage for the whole PE S region, as can be estimated from the 
spectrum of Figure 1 in lJo et aLl (1201 lb . The contribution of 
H2 fluorescent emission amounts to ~35% of the continuum 
intensity near the western edge of the lower part of arc A. The 
intensity levels of this region decreased from —2,000 CU to 
— 1,200 CU after subtraction of H2 emission. Two-photon ef- 
fects of hydrogen atom are clearly seen in the two strong Ha 
filamentary structures, arc A and arc B, whose centroids are 
depicted as thick black solid lines in FigurefTJa). 

As Figure [Tta) generally shows high FUV intensity in the 
regions where the color excess E{B-V) is large, we made a 
pixel-to-pixel scatter plot of the FUV intensity against dust 
extinction in Figure [Ttb) to check the correlation between 
them. A correlation is seen for E{B-V) < 0.14, especially for 
the FUV intensity less than 1500 CU, while FUV intensity 
becomes scattered and more or less independent of dust ex- 
tinction for E{B-V) above 0.14. In fact, it should be noted 
that the color excess value of —0.14 corresponds to the opti- 
cal depth of T — 1 at 1565 A. The result agrees well with the 
general notion that the FUV continuum and dust extinction 
have a correlation in an optically thin region. The data points 
of E(B-V) < 0.14 with its FUV intensity of less than 1500 CU 
correspond to the blue region in Figurellla): among these data 
points, those of higher FUV intensity correspond to the inner 
region of the OES shell while those of lower FUV intensi- 
ties correspond to the outer regions of the OES shell. Both of 
them seem to show a good linear relationship, implying more 
or less constant radiation fields. 

The data points of E(B-V) < 0.14 with its FUV intensity 
of above 1500 CU generally correspond to the regions of 
light blue to yellow around the boundary of the thick dust 
region. These regions show high FUV intensity due to the 
strong UV radiation fields produced by the nearby B type 
stars. The strong UV radiations from these bright B type 
stars also affect the optically thick region in the eastern part 
of the Orion Molecular Cloud Complex (henceforth OMCC): 
the data points of E{B-V) > 0.14 and the FUV intensity above 
2500 CU correspond to the regions of yellow to red in the 
eastern part of the Complex in Figure (TJ a). The data points 
of E(B-V) > 0.14 and FUV intensity < 2500 CU correspond 
to the regions of ligh t blue to green in the western part of the 
OMCC. As noted in iJo et all d2011l) . the western part of the 
Complex shows lower FUV intensity than its eastern coun- 
terpart because it contains fewer stars as interstellar radiation 
sources. Hence, we see a linear relationship between FUV in- 
tensity and dust extinction in the blue regions of Figure [Tta) 
where dust extinction is small while the scattered data points 
in Figure [Ttb) seem to be related to the fluctuations in the 
background UV radiation fields. 

3. SIMULATION MODEL AND RESULT 
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Fig. 2. — (a) FUV continuum map obtained from dust scattering simulation 
for a single dust slab, and (b) a pixel-to-pixel scatter plot of the FUV inten- 
sities for the observed and simulated maps in Figure [TJa) and Figure |2ja), 
respectively. 194 stai's selected as photon sources ai'e marked in (a): black 
symbols correspond to the stars closer than 200 pc, and red symbols corre- 
spond to the stars located farther than 200 pc. The symbol size is proportional 
to the 1565 A band flux of the TD-1 star catalog in log scale. 



As the observed FUV intensity generally follows the dust 
extinction level in Figure [Ua), with fluctuations which may 
be attributable to the local background UV radiation fields, 
most of the FUV diffuse emission in the OES region seems 
to originate from dust scattering of the radiation fields pro- 
duced by nearby stars. Hence, it might be possible to ob- 
tain useful information regarding the distribution and scat- 
tering properties of dust from the observed FUV continuum 
map by comparing it with simulations based on theoretical 
predictions. We have performed Monte Carlo simulations, 
in which photons originating from bright stars are multiply 
scattered by assumed distributions of dust. We adopted the 
following Henyey-Greenstein scattering phase function $((?) 
dHenvev & Greensteinlll94ll) for the dust scattering model, 
which is relatively simple and known to fit the observations 
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Fig, 3. — 90% confidence contours: (a) thickness vs. distance plotted for the albedo 0,43 and the g-factor 0,45, (b) g-factor vs, albedo plotted for the distance 
80 pc and the thickness 70 pc, (c) thickness vs, albedo plotted for the distance 80 pc and the g-factor 0,45. (d) distance vs. g-factor plotted for the thickness 70 
pc and the albedo 0.43, (e) distance vs. albedo plotted for the thickness 70 pc and the g-factor 0.45, and (f) thickness vs. g-factor plotted for the distance 80 pc 
and the albedo 0.43. Here, distance indicates the distance to the front face of the slab. 



reasonably well. 

^ 47r[l+^2_2^cos(6')]i-5 

where a and g are the albedo and phase function asymme- 
try factor of dust grains, respectively, and 6 is the scatter- 
ing angle measured from the direction of the incident pho- 
ton. The simulation cod e adopted the peeling-off method 
dYusef-Zadeh et all 119841) : while each photon experiences 
multiple scatterings in random directions as it propagates 
through the dust cloud, the code keeps track of the fraction of 
a photon in the direction toward the Sun from each scattering, 
which is summed up for the final intensity. The Mont e Carlo 
simulation code is described in more detail in ISeoiil (l2009h 
and Seon et al. (2012, to be submitted). 

We took a 400 pc x 400 pc x 800 pc rectangular box as 
a simulation domain, with 800 pc chosen to be the distance 
to the rear face from the Sun in view of the size of the OES 
region whose distance is known to be ^155 pc at the near- 
side and ^586 pc at the far-side. The central axis was chosen 
along the line of sight toward the center of OES, {a, S) = 
(60°, -5°), with the Sun placed at the center of the front face. 
The 400 pc x 400 pc square located at 750 pc from the Sun 
corresponds to 30° x 30° angular span in the sky, encom- 
passing the entire OES region of 30° x 30° of the present 
study. The bright stars of the TD-1 star catalog, whose 1565 
A band flux is greater than 10"'^ ergs A"' s"' cm"^, were 



selected as the point sources of stellar photons. The total 
number of stars is 3,338 for the whole simulation box, 194 
of which belong to the region of a 30° x 30° angular span 
of the present interest and are shown in Figure EJa). Incor- 
poration of the additional Hipparcos stars did not change the 
result much, though the Hipparcos catalog contains ^3,000 
more stars than the TD-1 catalog for the OES region: the total 
luminosity increase was only ^\%. The stellar luminosities 
were calculated with the distances adopted from the Hippar- 
cos star catalog, together with the extinction corrections made 
by comparing the observed (B-V) values and th e mo del (B-V) 
values d erived from the Castelli model (jPerryman et all 19971: 
ICastelli et al. 2003). We checked th e E(B-V) values fo r the 
19 OES stars listed in the catalog of lFriedemannI (119921) and 
found good agreement within 15% on average for sufficiently 
high E{B-V) as much as 0.02. The optical depth at 1565 A 
was derived from E{B-V) assuming Ri/ to be = 3.1, an av- 
erage value of the diffuse ISM, with the absorption coeffi- 
cients calculated for the carb onaceous-silicate grains model 
dWeingartner & Drainell200l]) . In fact, we estimated the aver- 
age Ry value for t he OES region using the values made for 
the 6 OES stars in Larson & Whittej (l2005h and the Ay val- 
ues provided by Neckel et al. (1980) for the 1 1 OES stars with 
E(B-V) > 0.05, and obtained 3.2, which was very close to the 
average value of the diffuse ISM. Furthermore, the extinction 
curves for five OES stars ( HD23 466, 24263, 26912, 30076, 
and 30870) in lPapai etalJdT99l are more or less consistent 
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with Rv =3.1. 

The most serious difficulty with the present simulation is 
related to the dust distribution since there is no report regard- 
ing the spatial distribution of dust for the OES region, other 
than the infrared survey results that give only indirect infor- 
mation on the integrated amount of dust along the sightlines. 
Hence, we assumed the dust cloud to be a slab of finite thick- 
ness for the sake of convenience. This idealized dust model 
may certainly be different from the real dust distribution, and 
the result of our simulations should be regarded as a represen- 
tative one, indicating the effective distance and thickness. In 
the present simulation, the E{B-V) values of the SFD dust map 
for given sightlines were converted into the optical depths at 
1565 A, which was regarded as a representative FUV wave- 
length of the FIMS observation, and then assigned to the grid 
points with equal distribution for an assumed cloud thickness. 
As we learned that most of the dust is located w ithin ^400 pc 
from t he Sun in the direction toward the OES dVergelv et al.l 
(120101) : Vergely 201 1 private communication), we took a sin- 
gle dust slab and placed it within 400 pc with the whole dust 
extinction from the SFD map distributed equally among the 
grids for an assumed thickness. We will further consider pos- 
sible improvements to our model in the Discussion section. 

Along with the distance and thickness of the dust slab, the 
albedo and phase function asymmetry factor, which are the 
optical parameters related to the scattering properties of dust, 
were varied. The ranges of simulation parameters were as 
follows: the distance to the front face varied from 10 pc to 
200 pc with 10 pc steps, the thickness varied from 10 pc to 
200 pc with 10 pc steps, the albedo varied from 0.35 to 0.50 
with 0.01 steps, and the g-factor varied from 0.20 to 0.70 with 
0.05 steps. For each simulation, 10^ photons in total were 
randomly generated from the point sources. Each simula- 
tion result was compared with the FIMS observation to obtain 
the best fitting parameters using the chi-square minimization 
method. 

Figure |2ja) shows the simulated FUV map for the best fit- 
ting parameters: 80 pc for the distance to the front face of 
the dust slab, 70 pc for the thickness of the slab, 0.43 for the 
albedo, and 0.45 for the g-factor. Figure |2b) is a scatter plot 
obtained from a pixel-to-pixel comparison between the FIMS 
map of Figure[Tla) and the simulation map of Figure|2ja). We 
have divided Figure |2j a) into nine sub-regions with a 10° x 
10° area and numbered the sub-regions from 1 to 9, so that 
the data points in Figure |2b) were color-coded according to 
the sub-regions. It can be seen in Figure |2b) that the simu- 
lated intensity generally agrees with the observed one. The re- 
duced chi-square value is ^ 5.44. A little lower simulated 
intensity than that observed below 1500 CU of FIMS inten- 
sity (blue) corresponds to the regions of arc B and the lower 
part of arc A, where two-photon emission, not considered in 
the simulation, is a significant contribution. If the difference 
is regarded as the contribution from two-photon emissions, 
the amount is estimated to be ^25-40% of the total observed 
co ntinuum , a lit tle lower than the previous value of '--^50-70% 
in lJo et ai] 1*201 1), which was obtained by converting the ob- 
served Ha intensity assuming that 1 Rayleigh of the Ha line 
corresponded to the two photon emission of 60 CU. 

Large fluctuations are seen between 1500 and 2500 CU of 
FIMS intensity (sky-blue), and the simulation result is consid- 
erably lower than the observed one above 2500 CU of FIMS 
intensity (orange). These discrepancies may originate from 
local variations of the simulation parameters such as the dis- 



tance to the front face and the thickness of the dust cloud. 
More consideration of these discrepancies will be given in the 
Discussion section. 

Figure |3] shows the contour maps for the 90 % confidence 
(Ax^ ^ 1.06) among the four variables, the albedo, the g- 
factor, the distance to the front face of the slab, and the thick- 
ness: the ranges are 0.43!I^'q4 for the albedo, OAS'^i for the 
g-factor, 80;!;5Q pc for the distance to the front face, and 70;!;gjj 
pc for the thickness. The 90 % confidence ranges for both the 
albedo and the g-factor include the t heoretical val ues 0.4 for 
the albedo and 0.65 for the g-factor ( lDrainell2003l) . It should 
be noted that the albedo is well constrained with less than 10% 
error range for the 90% confidence. The distance and thick- 
ness values indicate that the dust cloud is generally located in 
front of the OES, whose distance is known to be ^155 pc to 
the near-side. The wide ranges for the distance and thickness 
estimated here may actually indicate the large spatial dimen- 
sion of the cloud along the line of sight, as we will note in the 
Discussion section. 

4. DISCUSSION 

As mentioned previously, the assumption of a single slab 
as a model for the dust cloud in the OES is unrealistic, and 
the ranges of the estimated distance and thickness were rather 
broad. Here, we would like to improve the model by allow- 
ing local variations of the distance and thickness of the cloud. 
However, we fix the value of albedo to be 0.43, as we obtained 
with a single slab model. Being the most sensitive one among 
the four parameters that determine the model FUV intensity 
level, it was constrained very well. We also set the value of 
the phase function asymmetry factor to be 0.45. The range of 
g was found to be rather broad, implying that the model was 
not affected significantly by the change of g. Furthermore, 
Figure|3jd) and FigureQf), in which the 90 % contours of the 
distance and thickness were plotted against g, do not show any 
correlation with the value of g. With these fixed albedo and 
g-factor values, we carried out simulations for the sub-regions 
1 to 9 individually to obtain the best fitting parameters of the 
distance and thickness for each of the sub-regions. 

The best fitting values for the 9 sub-regions are shown in 
Table [T] and the resulting mosaic map constructed with these 
fitting parameters is depicted in Figure SJa), together with 
the pixel-to-pixel scatter plot of the FUV intensities in Fig- 
urelUb). As FigureHfb) shows, the simulation result is much 
improved now when compared with the previous result shown 
in Figure|2lb), with the new reduced chi-square value of x^ ^ 
3.19. The most remarkable change is seen in the sub-region 
7, of which the intensity increased. Noting that the thickness 
increased dramatically for this region, we believe the reason 
for the intensity increase in sub-region 7 is that more stellar 
sources reside in the cloud now and can contribute to the scat- 
tered continuum. We also note that the contribution from the 
two-photon effect, estimated as a difference of the FUV in- 
tensity between Figure [Tta) and Figure |4j a), has decreased to 
^15-30%, by ^10% from the previous result. 

While the best fitting parameters shown in Table [T]may still 
not be exact, they seem to provide meaningful insights into 
the morphology of the OES region. First, we note that the 
sub-regions 1, 2, 3, and 6, corresponding to the blue regions, 
have very thin (^10 pc) dust clouds located at '--^70-90 pc. As 
it is believed that inner X-ray cavities of hot gas confront the 
ambient dust clouds in the OES, the model result of thin dust 
layers in these sub-regions may be the front boundary of the 
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Fig. 4. — (a) A mosaic FUV continuum map obtained from dust scattering 
simulations, witli the distance and tliickness of tlie sub-regions varied individ- 
ually, and (b) a pixel-to-pixel scatter plot of the FUV intensities taken from 
the observed and simulated maps in Figure[Tla) and Figure|4la), respectively. 



OES shell encompassing the hot cavity. The clouds may have 
become thin as they were blown out by the stellar winds or su- 
pernova explosions that produced the ho t cavity. Sub-region 1 
includes th e MBM 20 molec ular cloud ( Magnani et al.lll985h . 
which Bur rows et al.l (Il993h placed at the in terface between 
the OES bubble and the Local bubble. While iJo et all (1201 Ih 
also estimated its location to be in front of the OES bubble 
because the cloud seemed to block the FUV emission from 
the OES, the present model does not constrain its distance to 
be nearer than the clouds associated with the OES in the sub- 
region, perhaps either because the number of pixels relevant to 
MBM 20 is too small to affect the simulation result or because 
there are no stellar sources that can discriminate its distance 
from the OES clouds. 

Next, we note that the clouds of the sub-regions 7, 8, and 9 
all have their centers at a distance '-^llO pc, as the distance in 
the table is actually the distance to the front face of the slab, 
confirming it to be a single body. Hence, this part of OMCC 
seems to be located at a distance ~110 pc from the Sun, with 



TABLE 1 

Best fitting parameters for each sub-region with fixed 
albedo 0.43 and g-factor 0.45 
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Fig. 5. — E{B-V) values derived from the dust models of the present 
study are compare d with those for the 19 OES stars listed in the catalog of 
IFriedemannI (1993) . 



its thickness ranging from 130 pc at the core (sub-region 7) 
to 50 pc at the boundary (sub-region 9). Sub-regions 4 and 
5 are located between the two distinct groups of dust clouds 
mentioned above and the simulation results show quite differ- 
ent thicknesses for the two sub-regions while their distances 
to the front face are similar. We believe such a discrepancy 
arose because two different types of structures, a thin shell 
surrounding the hot cavity and a thick dust layer as part of 
the OMCC, coexist in these sub-regions. In fact, when the 
simulation was carried out separately for the upper and lower 
parts of these sub-regions, their thicknesses were ^100-140 
pc for the upper parts and 10-30 pc for the lower parts for 
both of the sub-regions, in accordance with the environments 
of a thick dust region to the north and a thin shell to the south. 
In reality, the structure of the dust clouds would be more com- 
plex than we have assumed in the present study. However, we 
believe that our results provide an overall structure of the dust 
clouds toward the OES region. 

With the new distance estimation for each of the sub- 
regions, we would like to compare the E(B-V) values derived 
from the distances of the dust layers with the obse rved E(B-V) 
values for the OES stars listed in the catalog of IFriedemannI 
(119921) . While the derived E{B-V) values suffer from many 
error sources, as the present study modeled very simple dust 
layers and the observed E(B-V) values are based on the stellar 
models that involve uncertainties, we believe such a compari- 
son should give insight into the validity of our model. Figure 
|5] shows the result: the derived E{B-V) values plotted against 
the observed E(B-V) values. Although there are a couple of 
outliers, a positive correlation can clearly be seen, with a cor- 
relation factor of 0.75. However, the derived E(B-V) is gener- 
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ally higher than the observed E(B-V). While we may expect 
that real dust distributions are not as simple as in the present 
study, such a tendency may come from an over-estimation of 
the amount of dust placed in front of the stars: in fact, we put 
the whole amount of dust estimated from the SFD survey in 
the model layers, without leaving any for the background. 

Finally, it would be of great interest to know whether the 
characteristics of the dust in the OES region are comparable to 
those in previously studied regions, as the OES was suggested 
to have been created by an energetic supernova explosion with 
an additional energy source from the stellar winds produced 
by the stars in the I Ori association ( Reynolds & Ogdenri979h . 
and the origin of the OES might be reflected in the proper- 
ties of the dust. It has been shown that dust grains charac- 
terized by relative flat FUV extinction curves stand out with 
relati vely high FUV alb edo values, in the cases such as IC 
435 ( Calzetti et all [19951) and ScoOB2 (IGordon et al.|[T994h . 
while dust with steeper far-UV extinction curves, suggestive 
of a lar ger component o f small grains, exhibits lower albedo 
values (I Witt et al] 1 19971) . We fit each of the 9 sub-regions 
by treating the albedo as a free parameter to see if there is 
a difference between the thick dust region associated with the 
Orion Molecular Cloud Complex in the north and the cavity 
region where hot gas prevails. However, we could not find 
any systematic changes in the albedo among the sub-regions, 
with its values varying from 0.42 to 0.47. We also checked the 
Ry values to see if there is a systematic variation among the 
sub-regions using the data from Larson & Whittet (2005) and 
the Ay values provided by Ne ckel et al.i ( 1980), as mentioned 
in Section 3, but could not find a meaningful change between 
the northern molecular cloud region and the southern cavity 
region, although statistics may not be high enough as there 
are only two stars located in the cavity region. Hence, the 



dust characteristics obtained in the present study seem to be 
similar to those of the diffuse ISM, with no exceptional signa- 
tures in the albedo and the Ry values. We believe the reason 
is probably that most of the dust seen in the present study is 
located outside the OES. In fact, the structure of the dust layer 
estimated in the present study, with its thickness of 10 pc at a 
distance of 90 pc in the cavity regions, is in good agreement 
with this view of it as a thin layer between the OES and the 
Local Bubble. 

5. CONCLUSIONS 

We have performed Monte Carlo simulations for the dust 
scattering of FUV emissions in the OES region and compared 
the results with the diffuse emission map made from the re- 
cent FUV imaging spectrograph mission, FIMS. We were able 
to obtain the optical parameters of interstellar dust grains for 
the OES region: the average albedo is 0.43;|;Q g| and the phase 
function asymmetry factor is 0.45^q 2' agreement with pre- 
vious observational and theoretical estimations. Furthermore, 
the simulation results indicate that the dust clouds are located 
in front of the OES, with its distance of ^110 pc and thick- 
ness of ~'50-130 pc, while the hot X-ray cavity is bounded by 
a thin (^10 pc) dust shell toward the Sun. 
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